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We study the exact non-singular zero-surface tension solu- 
tions of the Saffman- Taylor problem [Mineev-Weinstein and 
Dawson, Phys. Rev. E 50, R24 (1994); Dawson and Mineev- 
Weinstein, Physica D 73, 373 (1994)] for all times. We show 
that all moving logarithmic singularities aj, (i) in the complex 
plane oj = e"'', where </> is the stream function, are repelled 
from the origin, attracted to the unit circle and eventually 
coalesce. This pole evolution describes essentially all the dy- 
namical features of viscous fingering in the Hele-Shaw cell 
observed by Saffman and Taylor [Proc. R. Soc. A 245, 
312 (1958)], namely tip-splitting, multi-finger competition, in- 
verse cascade, and subsequent formation of a single Saffman- 
Taylor finger. 



Saffman- Taylor experiment. In 1958 Saffman and Tay- 
lor [1] studied the displacement of oil by water in a nar- 
row gap between two parallel plates (the Hele-Shaw cell 
[2]). The three upper pictures in Fig.l show three consec- 
utive stages in their observations: initial, intermediate, 
and asymptotic. The small ripples in the left-most pic- 
ture are caused by an instability of the initially almost 
planar oil/water interface. We call this initial stage per- 
turbative since the deviation from a flat front is small, 
so the dynamics can be examined by perturbation the- 
ory [3]. Later, at the intermediate stage shown in the 
center picture, many fingers are formed along the mov- 
ing interface. They compete, some are screened and then 
stagnate, some continue to advance and fatten, while new 
ones continue to appear because of tip-splitting. We call 
this second (intermediate) stage non-perturhative, since 
the deviation of the front from a straight line is no longer 
small and thus it cannot be studied perturbatively. Even- 
tually, in the final asymptotic stage (the right picture) , a 
single uniformly advancing finger dominates after all oth- 
ers are suppressed in the competition described above. 
This single finger occupies one half of the Hele-Shaw cell 
width w and is described by the formula 



X = 2t + log (cos y) 



(1) 



where w = 27r in our scaled units; the longitudinal ex- 
tent, that is the x-direction, is taken to be infinitely 



large; and velocities of both fiuids are chosen to be 1 
when |a;| -^ cx). Throughout the last 40 years after 
[1] was published, this experiment has been reproduced 
many times, both in laboratories and by numerical cal- 
culations. (See [4] and references therein). The process 
above was found to be universal and independent of the 
specific liquids (oil and water in [1]) provided the fluids 
are immiscible and incompressible and a driving fluid is 
much less viscous than a driven one. 

A mathematical form,ulation in the absence of surface 
tension and in scaled units is 

V^p = (in the oil domain), 

p ~ —X if a; ^ -l-oo (oil pushed to the right), 

p = (at the oil/ water interface), (2) 

dnP = for y = ±7r (at the channel walls), 

Vn = —dnP (at the oil/ water interface), 

where V„ is the normal velocity of the interface, p is pres- 
sure, and dn is the normal derivative. From (2) follows 
the Laplacian growth equation (LGE) (see [5] and refer- 
ences therein), which can be written compactly as: 



Im(ztZ0) = 1 , 



(3) 



for the moving front z{t,(j)) = x{t,(j)) + iy{t,(j)), 
parametrized by G [0, 27r]. Here the bar denotes the 
complex conjugate; zt and z^ are partial derivatives; and 
the map (f) —^ z{t, (j)) is conformal for Im cj) <Q. 

Derivation of the LGE (3) from (2) follows. Because 
of the identity Vn — Im(ztZ/), where I is an arclength 
along the front, and of the Cauchy-Riemann relation, 
—dnP — di(j) for pressure p and stream function 0, 
where di is the tangential derivative, the last equation 
in the system (2) can be rewritten as 

Im(ztz/) = di4> 

which is equivalent to (3); {dicj) — —dnP ^ 0, since p = 
along the front and because of the maximum principle 
for harmonic functions). It then follows that the variable 
(j) in (3) is the stream function. The map (j) ~* z(t, (j)) 
is conformal for Im (/> < 0, since the complex velocity 
potential W{z) — —p{x, y) + i4){x, y) is analytic in the oil 
domain. 

Previous studies have focused primarily on the selec- 
tion of the asymptotic finger width A from the continuous 
family, fomid in [1]: 
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A)log(cos^) 



(4) 



where the finger moves to the right with velocity 1/A, and 
the finger's dimensionless width A, measured in units of 
the channel width, can be any number between and 1. 
(It is not difiicult to see that (4) is the traveling-wave 
solution of (3)). These selection studies focused on sur- 
face tension as the critical factor, which determines the 
finger width A [4,6]. (When A = 1/2, (4) becomes (1)). 
Some studies (mostly numerical and semi-analytical) also 
addressed the intermediate (non-perturbative) period of 
the process [7], and again, as in the case of width se- 
lection, all the main features observed during this stage 
were attributed to non-zero surface tension. We believe 
that this attitude was formed because almost all exact 
solutions known at that time [8][J for zero surface ten- 
sion exhibited finite-time singularities (cusps). The in- 
evitability of cusps was related to the ill-posedness of the 
problem in the absence of surface tension [10]. Recently 
however the selection of A = 1/2 was solved [11] in the ab- 
sence of surface tension using the exact time-dependent 
TV-logarithmic solutions of the LGE (3), which remain 
non-singular and analytic for all time (no cusps) [5] . 
These N-logarithmic solutions have the form [5]: 



N 



z{t, 4>) = T{t) + i^l(j} + ^ afc log(e*'^ - afe(i)) : 



(5) 



/c=l 



^N 



where a^ = const, |afe| < 1, /i = 1 — '}2k=i^k, and 
/i G [—1,1]. With some constraints on {a*:}, which are 
not specified here (see [12] for details), these solutions 
remain regular and do not exhibit finite-time singularities 
[5], unlike all other known time-dependent solutions of 
(3). The time dependence of afe(f) and T(i) is given by 

N 

/3fc = r - log flfe + ^ ai log(l - afeO/) = const (6) 
1=1 

N N 

and * + C' = 'r+^X!X! "'="' ^°S(1 - a,kai) (7) 

k=l 1=1 

where C and Pk {k — 1,2, ..., A^) are constants in time. 
The validity of (6) and (7) can be easily verified after 
substitution of (5) into (3) [13]. The solutions (5) are 
dense in the space of analytic curves [14], and all singu- 
larities of z^ are simple poles ak{t), while all other known 
solutions of (3) [8] have multiple poles. Thus (5) are the 
most natural and generic exact solutions of (3) with a 
finite number of moving singularities. 

The geometrical interpretation of (5) is of great help 
(see Figure 2): The constant of motion (3k — afclog2 is 
also the complex coordinate of the stagnation point at 



the vertex of the fc*'* "groove" with parallel walls. This 
groove originates during the front evolution; 'K\o.k\ is the 
width of this groove, and arg(afc) is the angle between the 
axis of symmetry of the groove and the horizontal axis. 
This interpretation of solutions (5) is in excellent agree- 
ment with all known low surface tension experiments; 
and we claim that all stages of this process are faithfully 
described by (5) (to be more exact - appropriate subset 
of (5) [15]). 

The goal of this paper is to show how the exact solution 
(5) describes essentially all stages of the Saffman- Taylor 
process, but especially formation of a single Saffman- 
Taylor finger. Incidentally, we also show how the zero- 
surface-tension solution (5) can reproduce the following 
dynamical phenomena observed in this process and at- 
tributed traditionally to surface tension: tip-splitting, 
side-branching, screening of retarded and coarsening of 
advanced fingers, and formation of a single finger after 
finger competition. 

To accomplish our goal we will study the dynamical 
system of the poles a^ (t) defined by (6-7) on the mathe- 
matical plane uj = e^"^ . As mentioned above, the dynam- 
ics of poles Ofe provides, through (5), a complete descrip- 
tion of these phenomena in the physical plane z — x + iy. 
In particular, the coalescence of all poles, which we prove 
below, explains why only a single finger survives in the 
long-time asymptotics. The pictures at the bottom of 
Fig. 1 show typical distributions of singularities Ok inside 
the unit circle associated with all three stages of this pro- 
cess. 




^Except for a few nonsingular examples [9], limited however 
to a specific symmetry, and hence non-generic. 



FIG. 1. Three consecutive stages of viscous fingering in 
the Hele-Shaw cell: initial (left), intermediate (center), and 
asymptotic (right). The physical plane z is shown in the up- 
per pictures, while the lower pictures depict a distribution of 
moving poles ak{t) in the unit circle \(jj\ < 1 on the mathe- 
matical plane u. The open circle indicates the repeller, tj = 0, 
while the solid circle indicates the attractor, cj = 1, of poles 
whose dynamics is given by (6-7). 

Initial stage: disordering begins, the repeller at uj — 0. 



For a front that is almost flat initially, all singularities are 
near the origin: \ak\ <C 1 (see the left of Fig. 1). Indeed, 
in this case the front is almost planar: z(0, 0) = t + i(f). 
Then the singularities a^ start to move away from the 
origin. One can easily see from (6) that, if \ak\ ^ 1, 
\o-k\ — e"^. This exponentially fast repulsion from the 
origin, holds for every pole a^, if |afe| <C 1, regardless of 
the positions of the other poles. It means that uj = is 
the repeller for the dynamical system of ak{t) specified by 
(6) and (7). As one can easily see, a^ = is an unstable 
equilibrium solution of the dynamical system (6-7). (The 
corresponding (3k equals infinity.) The repulsion from 
zero manifests itself in the front instability, which always 
occurs when a less viscous fluid pushes a more viscous 
one. Indeed, the closer Ofc comes to the unit circle, the 
more complex the interface becomes. 

Because of the exponentially fast instability, the sys- 
tem quickly enters into the non-perturbative stage: con- 
vex parts of the interface accelerate, while concave parts 
decelerate, retard, eventually stagnate and start to form 
grooves near their most retarded points (stagnation 
points). The tendency of moving fronts to form these 
grooves with almost parallel walls has been observed 
in numerous laboratory and computational experiments 
[4,6,7]. The dynamics of formation of these grooves co- 
incides with the dynamics of the solution (5), given by 
(6) and (7), near the points z{t,4>) = j3k — afelog2 and 
is fully compatible with the geometrical interpretation 
of (5) given above. It is unknown, though, why these 
grooves have parallel (or almost parallel) walls. 

Early intermediate stage: disordering increases. The 
formation of grooves near stagnation points explains such 
frequently observed dynamical features as tip-splitting 
and side-branching of the moving interface. Each such 
event makes the moving interface more complex: it in- 
creases the effective number of grooves by one, and conse- 
quently the total number of grooves grows in time. (One 
can equally say that the number of fingers grows in time, 
since a finger, as we traditionally call it, is an element of 
the advancing front between two adjacent grooves. Note 
that while grooves have parallel walls, fingers, as a rule, 
do not.) In short, disordering and complexity of the front 
increases during this stage (the center of Fig.l). 

As shown in [5], the only long time asymptotics for 
each ak{t) consistent with (6) is \ak\ — > 1. This was 
shown assuming that Reafe > for each k = 1,2,..., TV. 
This assumption was taken in [5] as a working criterion 
for the solutions (5) to be non-singular. Recent stud- 
ies showed that this criterion is too simple [12], but for 
many scenarios it works successfully. Since (6) and (7) 
require all singularities ak{t) to move toward the unit 
circle (see details in [5]), they eventually enter its nar- 
row vicinity and accumulate there. The a^'s continue 
to approach the unit circumference \lo\ — 1 exponen- 
tially slowly, namely as 1 — |afe| = c^exp (— r/Reafe), 
when Cfc > 0, Re a^ > 0, and when the poles Ofc are 



not too close to each other. Only when ak approaches 
the unit circle does the groove near the stagnation point 
/3fe — ak log 2 develop. In this sense the geometrical in- 
terpretation of the constants ak , Pk given above is appli- 
cable only in intermediate asymptotics, that is when \ak\ 
approaches one. 

Late intermediate stage: ordering begins. "The law 
of nonlinear superposition" for groove merging. It was 
indicated and discussed in [5,16] that when two poles ak 
are close to the unit circle, they are mutually attracted 
along the angular direction, cannot pass each other and 
must merge. This merging is complete only when r ^ cx), 
while for finite r, the distance between merged poles is, 
as one can see from (6), of the order of e"^"^, where 6 > 0. 
As two poles merge, say ai and a2, we equate them, so 
the total number of logarithmic terms in the sum in (5) 
decreases by one, and a new logarithmic term in this 
sum forms from the two terms labeled by 1 and 2. The 
coefficient of this new term clearly is ai + a2- In the 
physical plane this corresponds to the merging of two 
grooves into a single one. The width and the angle (to the 
horizontal axis) of the newly formed groove are 7r]ai-f 02] 
and arg(ai -1-012) respectively, while the widths and angles 
of two grooves before their interaction were Trjai], 7r]a2l 
and argai, argQ;2 respectively. An example is shown in 
Fig.2. 
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FIG. 2. Geometrical interpretation of the complex con- 
stants of motion au and Pk; fc = 1, ..., A''. 

This "law of non-linear superposition" of grooves gov- 
erns the interaction between these grooves and deter- 
mines all the dynamics of groove coalescence. Each in- 
teraction of this kind leaves one of interacting fingers 
behind, while the others continue to advance and fatten, 
thus "winning the competition". Two adjacent fingers 
do not interact only if they are strictly parallel. Later 
we will show that this parallel configuration is linearly 
unstable. As a result of this merging the total number of 
advancing fingers is decreased by one. (Or, what is the 
same, the total number of effective grooves decreases by 
one after each pair coalescence). In this fashion, the solu- 
tions (5) describe the phenomenon of finger competition 
(or groove coalescence). This ordering process eventu- 
ally dominates the disordering due to tip-splitting and 
side-branching (which is just the birth and development 
of new grooves). So, in the long-time asymptotics only a 



single finger is left to propagate (Fig.l upper- right). 

Physically, the merging of grooves is an inverse cas- 
cade, an ordering process that dominates after satura- 
tion. This self-organizing behavior of an unstable front 
is not exceptional: it occurs in several unstable nonlinear 
processes in confined geometries [17] . It is the walls of the 
box the system is contained in, which stabilize growth. 
This happens when the transverse correlation length, ini- 
tially small, becomes comparable with the lateral width 
of the box [18]. Then in all these processes a univer- 
sal pattern appears (different for different processes). In 
the dynamics described by (2), this pattern is the single 
finger, observed by Saffman and Taylor in [1]. 

The rest of this paper derives the coalescence of all 
singularities ak at a single point on the unit circle (for 
periodic boundary conditions at the walls) or at two op- 
posite points on the unit circle a/c — > ±1 (for no-flux 
boundary conditions at the walls). In the physical plane, 
this means that, for generic initial data in the class of 
non-singular solutions (5), the asymptotics i ^ oo is in- 
deed a single uniformly advancing finger in accordance 
with experimental observations. Then replacing every 
Ofc by 1 and t by Vt in (5) when t — > cxd, we obtain 

z{t, <fi) = Vt + fii^ + (1 - A*) log(e*'^ - 1) 

which coincides with (4) when V — 1/A and fi — 2X ~ 1 
and thus describes a single uniformly moving finger, as 
expected. (To select A = 1/2 within the class (5), it is 
sufficient to move the pole infinitesimally away from the 
unstable equilibrium that is zero, as was shown in [11].) 
The variable t is a scaled time. First we prove that 
r in (5-7) is but a scaled time, and it will be treated as 
such below. First, we note that dr/dt is always finite, 
since we consider only non-singular solutions (5). Then 
we see that always t > t. To prove this, we expand every 
logarithm in the RHS of (7) as an infinite power series. 
We obtain 



ak when t ^ cxd. First, we introduce new variables ak 



t^T- 



OO 

E 

m— 1 



lE^ 



^i«fc«rP 



(setting C = in (7).) Thus t > t. Then if t 
does T. It follows from (6) and (7) that 



t + C 



1 



1 " 

-^afelog(afc) 



(8) 



fc=i 



where C" is a constant. The only way to compensate 
the divergent positive term t in the RHS of (6) is to 
admit that |afc| -* 1 for each k = l,2,...,iV [19]. The 
only way to compensate the divergent positive term r 
in the RHS of (6) is to admit that \ak\ -^ 1 for each 
A; = 1,2, ..., A^ [19]. Therefore in the long-time limit the 
relation between t and t is linear. 

Proof that all poles coalesce in the long-time limit. We 
want to investigate the asymptotic behavior of the poles 



exp{-?7fc+i$fc} and at = a'k 



The condition \ak\ 



1 then simply means that rjk -^ 0+. In order to eliminate 
the divergent term log(l — |afc|^) from (6) we multiply (6) 
by Oik and then extract its imaginary part. The result is 

lm{akl3k) = a'^^k - a'^T 

N 



+ 22[{a'i'a'f. -a'^a'i)\og\l-akai\ 
+ (afe"/ + a/'a'fe) arg(l - ftfefl,)]. 



(9) 



(We neglected the small term, —a'lrjk). Here we have the 
divergent term — aj^'r. To compensate this divergence we 
need to assume that at least some poles coalesce. That 
is N poles Ofc are divided into N groups of merging poles 
ak so that |$fe — $/| — * for all members of each group. 
After summation of eqs. (9) over all poles within each of 
the N groups, we obtain 



N 

-I- > [(a" (51 



Q;„(5„)log|l -a„o„| 
+ (5™"^ + "1<5") arg(l - &„,an)] ■ 



/ J LV nm ^^m ni 



(10) 



Here Cm are constants, <I>m = arg a™ , a^ = Z]fc=i ^'k ' 
a'm — X]fe=i o^fc' '^m is the number of poles in the m-Vci 
group, and 1 < tti < A^. Also because of the partial 
merging, all poles within each group labeled by m are 
equated to cim with an exponentially high accuracy. Re- 
placing (with the same accuracy) \am\ by 1, and assuming 
no merging between different am's, which means 

where k and I are from different groups, we rewrite (10) 
as 



N 
+ E [("""m - ^'La'n) log I Sin( 

+ ("«"« + "nam) 7, h 



$rr, - ^r 



(11) 



where C^ are new constants. From this equation we 
conclude that 

(i) Since log|sin(i^^)| ^ -oo for |l>„ - $„[ ^ 
0, 27r, the poles cannot pass each other. 

(ii) After summation over all m = 1, ..., N ^ we see that 

Em "m^m = COnst. 

(iii) From (i) it follows that once < $m — $„ < 27r 
in the long-time asymptotics it will hold for all further 
times. 

(iv) From (ii) and (iii) it follows that \^rn\ — > oo is 
impossible. 



(v) From (iii) and (iv) it follows that there is no way 
to compensate the divergent term —a'l^T in (11), except 
to admit that a'^ = for each I = 1, ...,N. 

Finally we conclude that 



i=i 



0, l>™ == 0, $„ ^ $„ (for m^n) 



To find the asymptotic motion of poles within each group 
containing rim poles we differentiate the equations (7) 
keeping only leading (divergent) terms: 



^ + E 



ai 



1=1 



VI + Vk + i{'^k - ^i) 



The solution of these equations is 



Vk = 
^k 



rile-^"-'- 



^, 



{^l-'^'i)e 



-r/d^ 



provided that 



E 
fc=i 



(12) 



Instability of configurations with two or more parallel fin- 
gers. Therefore, the only case when merging will not 
be complete occurs when two or more groups of sums of 
Qffc's, that are cim's, are real. This corresponds to the case 
(see the geometrical interpretation of a's above) when all 
newly formed grooves characterized by Sm's are parallel 
to each other (and to the horizontal axis). In this case in- 
stead of a single finger in the asymptotics, we would have 
several parallel fingers with grooves between them having 
widths of TT \a.m\- But this configuration is unstable: In 
order to verify this, let us perturb slightly the constants 
a'l so that (12) is violated for every m — 1, . . . ,N, or 
we may merely add new logarithmic terms to the sum in 
(6). Once we destroy (even slightly) the parallelism of 
the newly formed grooves that are partial sums of afc's, 
then in accordance with the analysis above we have only 
one group of poles which satisfy (12). And this group 
contains all poles, so iV = 1. Hence we have total coa- 
lescence of all poles Ofc . In the physical plane this means 
merging into a single groove (or, what is the same, a front 
forms a single finger). 

The unique solution for t ^ oo for each pole is 



$fc = $0e-2*/(i-/^=) 
1 - /^^ > 



w = 1 in Fig.l). In the case of no- flux boundary condi- 
tionsQ which impose double periodicity^ with reflectional 
symmetry on this problem [11] (see also [20]), there are 
clearly two groups of poles mutually complex conjugate 
to each other. Because of this (mirror) symmetry, the 
attractors are real, and there are clearly two of them, 
cj = ±1, in this case, instead of one attractor as for peri- 
odic boundary conditions. 

Proof of total coalescence without any formulae. After 
the detailed analysis (10)-(13) leading to the necessity of 
total coalescence, we will show how to prove the main 
result of this paper without a single formula, by using 
only the geometrical interpretation of (5) and the "law of 
non-linear superposition" given above: Any two adjacent 
grooves whose axes of symmetry are not strictly parallel 
(which is clearly the generic case) will eventually coalesce. 
This will happen where their axes of symmetry cross. In 
accordance with Euclid, two non-parallel straight lines on 
the plane necessarily cross, and since each crossing means 
the coalescence of two grooves, eventually all grooves co- 
alesce. At first glance, this proof is valid only for those 
grooves whose axes of symmetry cross in advance of the 
moving front, and invalid when they cross behind it. But 
the proof still works because of the mirror symmetry im- 
posed by no- flux boundary conditions (see [11]). This 
implies that if the middle lines of two adjacent grooves 
cross behind the front, then one middle line and the mir- 
ror image of the other one necessarily cross in front of 
the interface, so again coalescence is unavoidable. 

Summary of the asymptotic stage. After a cascade of 
subsequent groove crossings and coalescence, new grooves 
will be formed, described by am-, which are sums of the 
afc's of the original grooves participating in coalescence. 
New grooves will coalesce by the same rules as the origi- 
nal ones, that is when their symmetry axes cross. Since 
the parallel configuration of two or several grooves is un- 
stable, eventually a single groove will be formed with a 
width of 7r(l — fi). Since this is a real number, this fi- 
nal single groove will be parallel to the Hcle-Shaw chan- 
nel's walls. Then the single finger (a moving front be- 
tween two grooves) will propagate along the Hele-Shaw 
cell in accordance with observations. The finger width is 
27r — 7r(l — fi) = 7r(l -I- fi), and the relative finger width 
measured in units of Hele-Shaw cell (chosen here to be 
27r) is A = 7r(l + /i)/(27r) = (1 + fi)/2. Regarding the 
selection of A = 1/2 from this family (that is, /i = 0) we 
will just repeat here the essence of [11], where two obser- 
vations were made: that the origin is a point of unstable 
equilibrium (the repeller) for singularities, and that all 
fingers with a width different from one half correspond 



(13) 



Here we used from (8) that (1 -I- fi)T = 2t for t ^ oo. 
Thus for periodic boundary conditions all poles Ofc coa- 
lesce at a single point on the unit circle (chosen to be 



^expressed in (2) with dnp = at the channel walls 
''The channel width becomes tt in order to preserve the so- 
lutions (5). 



to the case when one of the singularities Hes at the origin. 
Once this singularity is perturbed (released from zero), 
it moves with the rest of poles toward the attractor lying 
on the unit circle [11]. Then, once there are no singular- 
ities at zero, the asymptotic finger has a width A = 1/2, 
in accordance with experiments [1,4]. 

General summary. The solutions (5) indeed explain all 
stages of front dynamics in the Hele-Shaw cell in great 
detail. On the mathematical plane, the singularities Uk 
move away from the repeller at the origin, where most of 
them are initially located, toward the unit circle which 
is the intermediate attractor for the dynamical system 
(6,7). On the physical plane this is a process of disorder- 
ing which occurs in agreement with the geometrical in- 
terpretation, that is through the formation of grooves de- 
scribed by constants {a/c, /3fc}, until the saturation time 
marks the end of the disordering part of the process. 
Then the system enters into the ordering part, when the 
number of competing fingers eventually decreases and fi- 
nally a single finger survives the competition and per- 
sists to advance in the long-time asymptotic limit (Fig.l 
upper-right). On the mathematical plane it means that 
when the radial component of poles' velocities becomes 
small near the unit circle, the angular component brings 
pairs of singularities together, and they attract and co- 
alesce. Then groups of coalesced singularities move as 
units near the unit circle and eventually merge together 
at a single point on the unit circle (w = 1 in Fig.l), 
which is the only attractor in periodic boundary con- 
ditions. The no-flux boundary conditions imposing the 
reflectional symmetry make two attractors at ±1. 

Conclusion. We have shown that the front expressed 
by non-singular solutions (5), which describe Laplacian 
growth in the Hele-Shaw cell in the absence of sur- 
face tension, forms a single finger uniformly propagat- 
ing along the Hele-Shaw cell asymptotically in time. 
Essentially all the observed dynamical features of tip- 
splitting, finger competition, screening, coarsening, and 
finally forming a single finger are shown to be conse- 
quences of the solution (5) for a generic choice of initial 
parameters. 

What is unclear. Here we mention questions for further 
elucidation: We have studied the asymptotics far beyond 
the farthest stagnation point Pk, implying that for some 
reason perturbations of the front eventually stop. What 
physically happens is that during the formation of the 
single finger, the front gradually stabilizes and eventu- 
ally is no longer vulnerable to weak noise. Clearly, this 
stabilization cannot be justified without non-zero surface 
tension. It is remarkable that the zero surface tension 
solution (5) nevertheless works and explains all these 
phenomena including tip-splitting, side-branching, finger 
competition, screening and coarsening [5], formation of 
a single finger (the present paper), and the selection of 
its width of one half [11], in accordance with numerous 
observations. Why does it work without surface tension? 



And how does surface tension regulate population and 
location of stagnation points and width of grooves? The 
answers to these questions are not known today. 
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